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During the time since the Initiation of this grant, the major 
activity in support of this grant has been Dr. B. R. Barkstrom' s 
work on the Ad Hoc Committee on the Earth Radiation Budget Measure- 
ment Process, with emphasis on sensor calibration and measurement 
accuracy. Dr. Barkstrom' s contribution to this committee involved 
a significant fraction of his time since the inception of the com- 
mittee on November 7, 1978. The work included the preparation of 
summaries of the past work on the earth's radiation field that must 
be used in reducing observations of the radiation field. Preliminary 
portions of the committee's final report are included here as appen- 
dices, and represent the fo”mal portion of Dr. Barkstrom' s contribu- 
tion to this committee. 

Dr. Barkstrom also provided significant contributions to the 
research in the Atmospheric and Environmental Sciences Division by 
aiding with the implementation of the finite difference radiative 
transfer algorithm for development of models of the angular and 
spectral dependence of the earth's radiation field for use in data 
reduction of the Nimbus 6 scanner data. He has also been de/eloping 
a simple version of the optimization algoritlim. 
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REVIEW CF EARTH AND SOLAR RADIATION ESTIMATES 
Energy Balance 

The simplest possible model of t*e Earth would be a sphere of constant 
and uniform albedo, isolated in outer space, illuminated by the sun, and 
emitting its own thermal energy into the void. In a steady-state condi- 
tion, the enf rgy (in W) accepted from the sun (given by the left side 
of equation l) is balanced by the emitted thermal energy (given by the 
right side of equation 1); that is, 
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In this equation, there are four unknowns: R^, the radius of the Earth: 
a, the global albedo of the Earth; F^, the solar constant; and F^ , the 
emitted (infrared) flux. 

The (mean) radius of the Eerth R^ is known to be 6371 km from geodetic 

surveys. A generally accepted value for the extraterrestrial solar >rraj ia „ti I# 

—2 —2 — 
is 1360 Wm (or 2 cal cm ). The global albedo a is generally 

thought to be around 1/3, It would then follow from equation (l) that 

the (average) emitted thermal energy density i3 P’ = 227 Wm* . 

That an energy balance exists npnears to be a good assumption since 

life has existed on the surface of the Earth for a long period and only 

small changes in the energy balance would cause drastic changes in climates. 


However, estimates of global albedo and emitted thermal energy are 


uncertain. 


Winds and Climates 

Sven before a quantitative statement of global energy balance was 
possible in the closing decade of the nineteenth cent^jfy, the experiences 
of the European explorers led natural philosophers and scientists to 
speculate on the causes of the winds and climates of different regions 
of the Earth. Hadley in a famous paper in 1735 seems to have been the 
first to suggest an explanation of the large scale features that is 
still most commonly accepted. He suggested that solar heating ne’r the 
equator would force air upwards and that the cooling near the poles would 
force air downwards. By insisting upon conservation of angular momentum, 
the tradewinds could be explained. 

The impartant conclusion from this explanation is that atmospheric 
circulation is caused by imbalances between incident solar energy and 
the emitted thermal energy. The gross features of this imbalance are 
shown in figures 1 and 2 (from Lorentz, 1967). The latter figure illustrates 
how the northward or southward transport of energy depends on the latitudinal 
imbalance between the solar input and the terrestrial output. 

Atmospheric circulation is, of course, much more complex than north- 
ward or southward movement. The distribution of land and water is exceed- 
ingly uneven, and their radiation properties differ markedly. Thus the 
reflected and emitted energy vary, producing marked variations in the 
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energy balance of the Earth's surface, as illustrated in figure 3 (from 
Budyko, 1978). 

Although figure 3 shows an average amount of energy available to be 
converted to other forms, it does not show the expected variations (with 
time) in the energy driving both atmospheric ani oceanic circulation. Nor 
does it give any indication of energy available to evaporate water, which 
may fora clouds and precipitate elsewhere. 

Effect of Changes in Energy Balance 

« 

— p 

A decrease in the incident solar irradiance of 1% (about lU Wm ) 

would result in a 1.5°C decrease in the mean surface temperature of the 

Earth if cloudiness and earth albedo remain constant. If thi3 small decrease 

would become permanent, then the border of the polar ice cover would shift 

from 10° to 18° in latitude toward the equator, and the mean surface 

temperature of the Earth would decrease by 5°C. These changes correspond 

l 

roughly to those that occured during the quaternary glaciations [SMIC, 19Yl]» 

A change in the earth albedo of 0.01 (e.g. , from 0.30 to 0.29 or about 
lU Wm in the reflected radiant exitance) would lead to a 2.3°C change in 
the mean surface temperature due to redistribution of radiative sources 
and oinks [SMIC, 1971]. Problems resulting from the sensitivity of Earth 
conditions to small changes in the energy budget are treated in numerous 
publications, including the study Inadvertant Climate Modification 
(SMIC, 1971) from which these illustrations ore taken. 

These types of concerns have led to the summary given in Table 1-1 
as taken from the report Toward an Internationally Coordinated Earth Radiation 


3 


Budget Satellite Observing System: Scientific Uses and Systems Considera- 


tions (1978). (How about a short summary of Table 1-1 rel^vynt to ERBE ? 

I 'mini lit ") 

A Short History of Solar Radiation Estimates 
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A Short History of Earth Albedo and Radiation Estimates 


Before the advent of satellites, there were two major methods of 

F. yf* f 

estimating the Earth albedo. summarizes some of these estimates 

as well as those obtained from satellite measurements. Most material of 

the following discussion is taken from Levine [1967]. 

(1) In the meteorological method, the reflected solar lrradlance 
ii 

wasr assumed to be made up of three components: ground reflection, 

atmospheric reflection, and cloud reflection. By developing estimates 
of the fraction of the earth covered by clouds and their reflectivity, 
the global albedo could be found. 

That is, the reflected radiation is assumed to be the sura of three 
components: 

a F - [f a + (1-f) (a + a )] F 
e c e a © 


where u is the earth global albedo, F is the incident solur irradiance 
13 b9 

(about thOB Wm~), f is the fraction of the earth covered by clouds, 

% 

a is the cloud albedo, a is the albedo of the earth, and a is the albedo 
c e ’a 

of the atmosphere. The most widely accepted estimates were given by London 

in 1957: for n year's average, f z 0.51, » c 5 0.47, a p » 0.139* and 

a ss 0.086. Thus, a ~ 0.35, and 2/3 of the reflected solar radiation is 
e 

due to cloud3. 


ii t» A 

(2) Another method of estimating the Earth albedo by measurj[n?: the 
eorthsh'.ne reflected from the moon. This . J-y ;* — ■ « uj procedure has been used 
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by Danjon (1936, 195 1 *) and Bakos (I 96 U). Generally, the earthahine 

rtvJL 

measurements produce a larger albedo than that suggested by other^ mea- 
surements. Furthermore, the earth 3 hlne reaaurements suggest that the 
albedo is highest in the Northern hemisphere winter, whereas meteorological 
methods suggest that it is lowest then. 

The method has been criticized on a number of grounds, including (l) 
that the observations so far do not include large areas of the earth (since 
the observations used come mainly from France and thereby exclude the 
Pacific), (2) that the angular dependence of reflection upon time of year 
ha3 not been properly accounted for, and (3) that the measurement excludes 
the near infrared region of the spectrum. 

Clearly, the major diffico'ty in developing accurate estimates of 
the reflected component of the earth radiation budget is due to the 
uncertainties in the area of the earth's surface covered by cloud and in 
the reflectivity of the clouds themselves. Tor example, Robinson (1956) 
has estimated the reflectivity of an overcast sky to bo O. 65 , while 
Neiburger (19*<9) derived a value of 0.75, and Roach's aircraft measurement 
(1961) gave a value of O. 85 . This spread is characteristic of both ob- 
servations and theoretical estimates. 

(3) Discussion of spacecraft measurement results up to Nimbus 6. 
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EARTH RADIATION MODELS 


•A 

INTRODUCTION 

\ 

In its simplest form, the ERB experiments are intended to measure 

a simple quantity - the power per unit over input or exiting from the top 

of the atmosphere over time scales long enough for the power to have an 

Impact on the weather or climate. This simple view of the problem suffers 

from the fact that the top of the earth's atmosphere is being treated as 

though it were a sheet metal cover to the atmosphere. However, provided 

that the top of the atmosphere is takei. to be indicative of an altitude 

reasonably higher than the levels at which ozone absorption is important, 

the top of the atmosphere appears to be a well defined quantity. It 

should be noted that because a sphere with radius 6371 + 50 km has 
6 2 

3.2 x 10 km more area than a sphere of radius 6371 + 30 km, a constant 

power per unit area going through the smaller sphere will result in a 0.6% 

-2 

smaller power per unit area going through the larger. Thus, if 240 wm 

-2 

are going out through the smaller sphere, then 238.5 wm are going out 
through a sphere with the top of the atmosphere at 50 km. 

Aside from the difficulty in defining the top of the atmosphere, the 
measurement of the reflected or emitted radiant exitance suffers from a 
more fundamental difficulty: the radiation emerges not in one direction 

but *n many. As a result, it is necessary to reduce the number of degrees 
of freedom by developing earth radiation models of spectral, angular, 
spatial, and temporal variations. 

In quantitative work, the basic variation in angle and wavelength is 
described by the radiance, L^. If we consider a given (theoretically 

• «. . • • 
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miniscule) area on our (now well defined) top of the atmosphere, then 
the power carried away by light travelling within a given set of direc- 
tions and also within a narrow wavelength band is given by 

dE - L dA d.Q cosO dA, 
as i'lustrated by Figure 1. 

The radiant exitance is then found by adding together the contribu- 
tions to the power leaving dA and dividing by the area: 

■ /dn cos0 . 

Because represents a sort of average (or moment) of , it is con- 
venient for our purposes to define a function 

R(£) = L/M 

which describes the angular variation of the radiation field, but not its 
strength. For the shortwave reflected radiation, R is (to within a constant) 
known as the bidirectional reflectance. For the longwave radiation R is 
known as the limb darkening function. One of the major tasks of providing 
earth radiation models is to provide reasonable estimates of R. 

Unfortunately, there appears to be no hope of directly measuring M at 
the top of the atmosphere because we would have to integrate over all direc- 
tions at every point on this "surface". Because of the fact that radiance 
is constant along a path unless the atmosphere gets in the way, it can be 
shown through Gauss' Theorem that if we measure /-dAM on some surface surround- 
ing the TOA, we can simply take this power and reduce it to the TOA by taking 


the radio of the areas. In this sense, the problem reduces to that of 
simply sampling M on the outer sphere well enough so that we can be sure 
we have measured the total power accurately. This Is probably the basis 
for some of the early concepts of earth radiation budget measurement, as 
suggested by Figure 2. 

The trouble is that as our understanding of the problem has Increased, 
so has our desire for more detailed knowledge. We now hope to produce from 
measurements of the radiant exltance on the integrating sphere values of 
the radiant exitance as a function of position on the TOA. As we have 
already stated, our interest lies in pole-equator radiation balances and 
in longitudinal varlaticns in M. 

When we wish to quantitatively relate the measurement to M, we must 
return to a geometry similar to that shown in Figure 1, in which we relate 
the power emitted by an area dA which is received by an area dA' on the 
detector: 

dA'cosO* n 

dE ■ L.dX — cos0 dA . 

A P Z 

This relationship constitutes the fundamental "law" in analyzing the re- 
lationship between the satellite signal and the quantity we desire. 

Not all areas on the earth are dark, of course, so the total power 
falling on the detector involves an integration over the area of the top 
of the earth's atmosphere visible to the detector. The detector converts 
this incident power into a signal with a sensitivity S,(0',$'), so that the 
part of the output signal coning directly from the earth is 

ji - Id A /dX /JA' S. [— ) R.M . 

visible detector 
earth 

. * are-% 


area 


Although it may not be apparent from the form of the Integration, the 
integral over dA is equivalent to an integration over the f ield-of-view. 
Although a treatment of other Influences and errors is given elsewhere 
in this report, this equation should be considered to be the fundamental 
equation relating a given measurement m to the quantity desired M. It 
applies to both narrow and wide f ield-of-view instruments. Of course, 
for a NFOV scanner channel, the integration over the earth's surface be- 
comes small enough that we usuallv expect to be able to assume both R and 
M are constant across it, so that they can be removed from the integral as 
constants. 

Although the product R^M could as easily have been written (the 
radiance), we have used directional-spectral function R^ because of the 
large number of quantities that would otherwise be required. As we have 
tried to make abundantly clear, our real interest is in M as it varies 
over the top of the atmosphere. This variation may readily be described 
by assuming that the earth's TOA is divided into a number of tiles of 
roughly equal area, e.g. 250 km x 250 km (or 2 1/2° lat x 2 1/2° long, 
at the equator). Table 1 shows the number of tiles required for various 
resolutions, as well as the number of L values that would be required to 
give the same resolution in direction of propagation as is required in 
space 
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Resolution 

10*xl0* 

(1000km x 1000km) 

5*x5* 

(500km x 500km) 

2 1/2* x 2 1/2* 
(250km x 250km) 

No. of M 
Values 

100 

400 

1600 

No. of L • 
values for 
identical 
resolution 
in M 

5000 

80000 

1280000 


TABLE 1 


Although the construct ion of M would be possible from observations of L, 
the number of quantities to be determined is clearly enoi.gh to clog the 
computer when 2 1/2° x 2 1/2* resolution is desired, even if enough values 
can be sampled to perform the integration empirically. 

Physics of the Angular and Spectral Variations 

It is clear that the quantity R is important in relating the measure- 
ments to the radiant exitance that vj desire. Our success in reducing the 
data is clearly going to be related to our ability to determine it and to 
characterize it by as few numbers as possible. Tn a broad view, R contains 
four possible types of variations that depend upon tlte nature of the under- 
lying surface, as outlined in Table 2. 

In order to explain briefly the data summarized in Table 2, it will be 
necessary to comment briefly on the nature of the interactions between the 
earth, its atmosphere, and the reflected or emitted radiation. 

In the shortwave portion of the spectrum, from about 2G0nm(.2wm) to 
perhaps 3000-5000nm (3-5wm), the light impinging on the detector is reflected 
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sunlight. The solar spectrum has been studied Intensively, and numerical 
values of the monochromatic lrradlance from the sun are available from 
several sources (e.g. Thekekera, 1971; Labs and Leckel). Since most ob- 
servations of the solar spectrum have been made at the ground, the published 
values tend to cluster together In the visual spectral range (roughly 
0.35'rJm to 0.7um) and to be somewhat more uncertain In the UV (owing to 
ozone absorption) and In the near IR (owing to water vapor absorption). 

Over clear scenes, the sunlight is reflected from the ground and 
scattered by the atmosphere. Ozone absorbs almost all incident sunlight 
with wavelengths shorter than about 0.35 urn. In unpolluted air, the light 
scattered from the atmosphere has been deflected by Rayleigh scattering 
from molecules, so that the reflected light is Inversely proportional to 
the fourth power of the wavelength. In addition, because the scattering 
occurs In an optically thin atmosphere, the reflected light is strongly 
limp brightened. To a fair approximation, the radiation emerging has a 
radiance distribution that is approximately given by 

T X 4 

L, “ M. (n— 2 — jr (~) + a e ] . 

X Xo 4ttcos0 X surf 

Atmospheric haze will tend to decrease both the wavelength dependence and 
angular variation in the reflected radiance. The major complicating factor 
is the surface reflection, which constitutes about one-half of the reflected 
solar radiation. Over the ocean, the reflection may assume a nearly 
specular character, whose half-width depends directly on wind speed 
(Kattawar, Plass, and Guinn). 


Completely overcast scenes or snow, on the other hand, have reflection 
that Is dominated by multiple scattering from nearly nonabsorbing particles. 

As a result, except for absorption by water vapor in the near IR bands, the 
reflection is not strongly dependent on wavelength. The angular variation 
in this situation is more difficult to characterize, although it may b 
treated theoretically - albeit with considerable caution owing to the strong- 
ly forward peaked scattering functions of large particles. General refer- 
ences include Pattrldge and Platt, 1976; Travis and Hansen, 1975; Irvine. 
Roughly speaking, when the sun is overhead, the reflected radiation should 
be reasonably close to Lambertian. As the sun moves towards the horizon, 
the reflected radiation tends to become more strongly forward peaked since 
the light is reflected a':er one or two scatterings and still carries much 
information about its interactions with the water droplets or ice crystals 
making up the scattering medium. 

Between clear conditions and completely overcast or snow lies a 
dif ficult-to-characterize area of broken or scattered cloudiness (or patchy 
snow). Although the reflected radiation might hopefully be characterized as 
composite between clear and cloudy, the situation is likely to be considerably 
more complicated. Theoretical studies of radiative transfer in the atmosphere 
have only recently advanced enough to consider the influence of cloud geometry 
upon reflection (Cox and McKee, 1975; Davies, Arduini and Barkstrom; and 
others). It is known that after sunlight enters clouds of finite horizontal 
extent it may leak from the sides. As a result, the larger clouds tend to 
be more reflective (Reynolds and McKee, 1978) a fact that has boen verified 
experimentally. In addition shadows are introduced that can be seen from 


some directions and not others. Because cloud cover, cloud size, and 
geometry change markedly with time of day, weather conditions, and oeason, 
the reflection from clouds constitutes the greatest uncertainty In charac- 
terizing the reflection function. 

In the longwave portion of the spectrum, from about 5 Mm to beyond 
100 pm, the light Impinging on the detector Is light emitted by the earth, 
the gases In the atmosphere, or clouds. The basic features of the spectrum 
and the angular distribution are easy to understand If we recall that Is any 
wavelength band, hot bodies emit more energy than cold bodies. Since the 
temperature of the atmosphere declines as we go from the ground to the 
tropopause and becomes reasonably constant thereafter, in transparent 
regions of the spectrum, emitted power is high. When the atmosphere becomes 
opaque In a given spectral band radiation emitted by the ground does not 
reach the detector outer space. Radiation that does emerge comes from high 
in the atmosphere and has a lower power than an equivalent high temperature 
source. [e.g. Kunde, et. al., 1960] 

Because a light ray propagating obliquely through the atmosphere en- 
counters more material tnan does radiation travelling vertically, radiation 
that arrives at a satellite detector after travelling obliquely through the 
atmosphere started higher than did normally propagating radiation. As a 
result, less radiation comes from the limb than from the nadii-the atmosphere 
appears to be neatly universally limb-darkened. 

The major effect of clouds is to raise the lower boundary of the at- 
mosphere, or to make the "ground" seem colder. This fact is routinely used 
to detect clouds in temperature rctrcivals from satellites. Thus, the 


spectral variation in the longwave region is reduced, but the angular 
variation is not expected to be strongly affected. 

In short, the physics of reflection and emission has long been stu- 
died. The major uncertainties in our understanding lie in the areas 
covered by broken clouds. Unfortunately, it is difficult to overcome this 
lack of understanding by making observations as we shall see. 


